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ABSTRACT 


A  theoretical  analysis  is  presented  for  the  problem  of  the  trans¬ 
ient  growth  of  the  ice  shell  that  forms  when  water,  initially  at  the 
freezing  temperature,  flows  through  a  pipe  which  is  convectively  cooled 
to  a  temperature  below  the  freezing  temperature.  A  solution  is  sought 
as  a  perturbation  expansion  about  the  quasi-steady  state  solution. 

The  first  order  closed  form  solutions  are  found  which  reveal  the  in¬ 
fluence  of  external  convection  and  of  sensible  heat  in  the  ice.  The 
influence  of  non-uniform  convection  on  the  interface  shape  is  displayed 
graphically  for  two  representative  cases  and  the  resulting  effect  on 
the  pressure  gradient  is  calculated  by  means  of  an  auxiliary  analysis. 

Experimental  results  are  presented  for  the  case  of  superheated 
water  flowing  through  a  uniformly  cooled  pipe  which  is  inclined.  The 
results  for  the  limiting  case  of  zero  superheat  agree  well  with  pre¬ 
viously  obtained  results. 
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CHAPTER  I 


INTRODUCTION 


The  phenomena  of  solidification  and  melting  are  seen  in  many  pro¬ 
cesses  ranging  from  the  natural  formation  of  ice  on  lakes  and  rivers 
to  the  casting  of  metals.  An  understanding  of  these  phenomena  is  im¬ 
portant  to  the  engineer,  especially  in  Canada,  where  development  of 
northern  regions  is  assuming  major  proportions. 

One  problem  particularly  worthy  of  attention  is  the  transient 
solidification  of  a  liquid  flowing  through  a  tube,  the  wall  temperature 
of  which  is  below  the  solidification  temperature  of  the  liquid.  As 
the  liquid  undergoes  a  change  to  the  solid  phase,  thermal  energy  is 
released  at  the  moving  liquid-solid  interface  and  must  be  removed  by 
conduction  through  the  solid  phase.  In  addition,  some  heat  may  be 
removed  from  the  flowing  liquid. 

This  type  of  problem  arises  in  many  common  situations:  for  exam¬ 
ple,  a  pipeline  buried  in  a  region  of  permafrost,  or  a  water  main  ex¬ 
posed  to  a  "cold"  meteorological  environment.  A  less  common  situation 
would  be  that  of  hydraulic  lines  in  a  space  craft  which  may  be  sub¬ 
jected  to  extreme  environmental  conditions. 

As  early  as  1916  Brush  [1]  discussed  the  principles  governing  the 
freezing  of  water  in  mains.  Although  he  presented  no  theoretical  or 
experimental  results  he  suggested  that  the  formation  of  ice  in  water 
mains  is  dependent  upon  the  temperature  of  the  water  and  its  flow  rate 
as  well  as  on  the  wall  temperature. 

In  1939  Pekeris  and  Slichter  [2]  considered  ice  formation  on  the 
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outside  of  a  long  pipe.  Their  work  employed  a  series  expansion  in 
which  the  zeroth  order  solution  corresponds  to  the  complete  neglect 
of  sensible  heat  in  the  ice.  Higher  order  terms  incorporated  the 
sensible  heat  which  was  assumed  to  be  a  small  fraction  of  the  latent 
heat. 

London  and  Seban  [3]  developed  a  general  approximate  method  for 
the  solution  of  one  dimensional  freezing  problems  for  several  geomet¬ 
ries  including  freezing  in  cylinders.  Neglecting  sensible  heat  effects 
in  the  ice,  they  used  thermal  circuitry  to  develop  closed-form  solutions 
for  both  first  and  third  kind  boundary  conditions. 

Kreith  and  Romie  [4]  presented  solutions  for  the  case  where  the 
interface  velocity  remains  constant.  They  considered  convective  cool¬ 
ing  for  cylinders  and  spheres  where  the  liquid  is  initially  at  the 
fusion  temperature. 

In  1962  Allen  and  Severn  [5]  used  relaxation  methods  for  the  case 
of  two  dimensional  ice  formation  in  planar  systems  and  one  dimensional 
freezing  in  cylinders.  They  considered  first  kind  boundary  conditions 
with  the  liquid  initially  at  the  freezing  temperature. 

In  the  same  year  Poots  [6]  applied  integral  methods,  similar  to 
those  employed  in  solution  of  boundary  layer  equations,  to  problems  of 
solidification.  The  results  for  the  circular  cylinder  were  compared 
to  those  of  Allen  and  Severn  with  satisfactory  agreement. 

Zerkle  and  Sunderland  [7]  considered  the  problem  of  superheated 
water  flowing  in  a  horizontal  circular  tube.  Their  theoretical  analy¬ 
sis  considered  steady  state  conditions  and  neglected  the  effects  of 
natural  convection.  Their  experimental  results  showed  considerable 
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deviation  from  the  theoretical  predictions  due  to  the  effect  of  natural 
convection  in  the  experimental  system.  The  work  of  Des  Ruisseaux  and 
Zerkle  [8]  is  an  extension  of  that  of  Zerkle  and  Sunderland  and  is 
concerned  with  predicting  conditions  under  which  hydraulic  systems  will 
freeze  shut. 

Gort  [9]  studied  the  problem  of  freezing  inside  a  long  vertical 
cylinder  with  and  without  a  forced  internal  flow.  He  presented  approxi¬ 
mate  solutions  for  arbitrary  wall  temperatures  for  the  case  of  a  liquid 
at  the  fusion  temperature.  Experimental  results  were  presented  for  the 
case  of  superheated  water  flowing  in  a  vertical  tube. 

O’zisik  and  Mulligan  [10]  presented  an  analytical  investigation  of 
the  freezing  of  a  liquid  flowing  in  a  circular  tube.  Constant  wall 
temperature,  slug  flow  velocity  profile,  and  quasi-steady  state  heat 
conduction  in  the  solid  phase  were  assumed.  The  temperature  distribu¬ 
tions  of  both  phases  were  found  and  these  solutions  were  coupled  through 
the  heat  balance  at  the  interface  to  give  a  single  ordinary  differential 
equation  for  the  interface  location.  Thus  interface  profiles  were  found 
as  a  function  of  time  and  axial  locations. 

Freeborn  [11]  investigated  the  problem  of  superheated  water  flowing 
through  a  convectively  cooled  vertical  cylinder.  The  analysis  was  broken 
into  two  parts;  an  ice  free  zone  near  the  entry  where  the  superheated 
water  undergoes  some  cooling,  and  a  zone  in  which  an  ice  shell  forms  in 
the  cylinder.  The  method  required  the  solution  of  the  temperature  dis¬ 
tribution  in  the  liquid  and  solid  phases  with  the  energy  balance  at  the 
interface  providing  a  differential  equation  governing  the  interface 
location.  Experimental  results  were  presented  with  satisfactory 
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agreement  between  theory  and  experiment. 

The  purpose  of  this  thesis  is  to  study  the  problem  of  asymmetric 
freezing  in  convecti vely-cooled  circular  tubes  under  transient  condi¬ 
tions. 

A  theoretical  analysis  is  presented  for  fully-developed  ice  forma¬ 
tion  in  which  the  effect  of  sensible  heat  in  the  ice  is  incorporated 
as  well  as  the  effect  of  non-uniform  external  convection.  The  effect 
of  asymmetric  ice  growth  on  the  internal  flow  is  calculated  by  an 
auxiliary  analysis. 

Experimental  results  are  presented  for  a  different  asymmetric  sit¬ 
uation  in  which  superheated  water  flows  through  inclined,  convectively- 
cooled  cylinders.  Axial  as  well  as  azimuthal  profiles  of  the  inter¬ 
face  are  presented. 


PART  I 


THEORETICAL  ANALYSIS 


CHAPTER  II 


ICE  GROWTH  ANALYSIS 


2.1  Governing  Equations 

Consider  water,  at  the  fusion  temperature,  flowing  through  a  very 
long  circular  tube  which  is  suddenly  immersed  in  a  coolant  (see  fig. 
2.1).  The  differential  equation  governing  conduction  in  the  ice  shell 
that  forms  in  the  tube  is  given  by 

32T  1  3T  1  32T  32T  _  1  3T  0  ,  , 

3R2  R3R  R2  3G2  3Z2  '  K3t  ’ 

where  T(Z,R,0,t)  is  subject  to  the  conditions 

KZ.R.0,0)  =  Tf  , 

T(Z.A,0,t)  =  Tf  , 

and  ki  |^(Z,R1>0,t)  =  -h(T(Z,R^ ,0,t)  -  Tro)  . 

We  introduce  the  normalized  non-dimensional  variables 


r  =  R/R-j  , 

0  =  0/  2tt  , 


z  =  Z/Zs  . 
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FIG.  2.1:  COORDINATE  SYSTEM 
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T  =  t/t  , 


and 


=  A/R1  , 


where  Zg  and  t  are  length  and  time  scales,  respectively.  The  charac¬ 
teristic  length  scale  for  the  problem  is  chosen  to  be  the  length  of  the 
tube.  A  characteristic  time  scale  may  be  found  in  terms  of  the  other 
characteristic  quantities  by  estimating  the  time  required  to  remove  the 
heat  released  in  completely  solidifying  the  water.  The  time  scale  thus 
calculated  for  the  mixed  or  third  kind  boundary  condition  (finite  h)  is 
not  the  same  as  that  found  for  the  Dirichlet  or  first  kind  boundary 
condition.  Since  the  Dirichlet  boundary  condition  is  a  limiting  case 
(h  ->  oo )  of  the  mixed  boundary  condition,  a  time  scale  was  chosen  that 
is  characteristic  of  a  mixed  boundary  condition  for  small  h  but  for  very 
large  h  is  characteristic  of  a  Dirichlet  boundary  condition.  Thus  we 
take 

t  : 


PiLR1  1  +  Bi. 


k. (Tp  -  T  ) 

i '  f 


Bi 


Substituting  the  normalized  variables  in  equation  2.1-1  gives 


dr  r  r  4  it  r  36^  l/  3z 


3t 


where  a  = 


Ste^  Bi^ 
"1  +  BTT 


This  simplifies  by  virtue  of  the  fact  that 


4tt  »  1  and,  for  a  very  long  tube,  Z$  »  R-j ;  and  therefore,  provided 
r  is  not  too  close  to  zero  it  follows  that 


Ol  +  1  H.  =  rv  M. 
9r2  r  9r  a  3x 


2.1-2 


9 


offers  a  good  approximation  to  the  description  of  conduction  in  the  ice. 
The  normalized  initial  and  boundary  conditions  are  then: 

cj)(r,  0,  a)  =  0 

4> (a ,  t,  a)  =  0  2.1-3 

and  (1  ,  x,  a)  =  -Bi.  [1  +  <j>(l,  x,  a)]  . 


For  this  radially  symmetric,  fully  developed,  situation  the  rate 

at  which  heat  is  released  per  unit  length  by  the  moving  ice-water  in- 

dA 

terface  is  equal  to  -Lp^  2tt  A  .  Assuming  the  water  remains  at  the 
fusion  temperature,  there  will  be  no  heat  transferred  to  or  from  the 
water.  Therefore,  the  heat  liberated  at  the  interface  must  be  conducted 
through  the  ice.  Thus,  radial  interface  motion  is  governed  by 

q  I  k  —  (A  t) 

pi  L  dt  i  3R  ^  j  ’ 


The  normalized  form  of  this  heat  balance  equation  is 


da 

di 


the  solution  of  which  must 


1  +  Bii  a* 

~TT.  8r  T* 


satisfy  the  condition 


9 


2.1-4 


a(0,  a)  =  1 


2.1-5 


if  the  system  is  to  be  initially  ice-free. 
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2.2  Solution  of  Equations 

In  a  great  many  situations  Ste^  «  1  and  since  a  is  always  less 
than  Ste.j  it  seems  reasonable  to  seek  a  solution  in  the  form  of  a  per¬ 
turbation  expansion  in  powers  of  a.  Thus 

00 

4>(r,  t,  a)  =  A  (r,  t)  +  £  a11  <f>  (r,  t)  .  2.2-1 

0  n=l  n 


Likewise,  a  solution  to  equation  2.1-4  may  be  sought  in  the  same  form: 
that  is 

00 

a(x,  o)  =  a  (t)  +  l  an  a  (t)  .  2.2-2 

0  n=l  n 


Substituting  2.2-1  into  2.1-2  and  grouping  like  powers  of  a  we 
obtain  an  infinite  set  of  differential  equations  for  the  temperature. 
The  first  two  of  these  are 


2.2-3 


and 


32<J>1  •,  3<J>-,  _  34>0 

3r2  +  r  "5F  "  57" 


2.2-4 


which,  to  satisfy  equations  2.1-3,  are  subject  to  the  boundary  and  ini- 
tial  conditions 

<J>0(r,  0)  =  0 

Vv0  =  0 

-gp-  (1  •  T)  =  -Bi.  [1  +  4>0(1 ,  T)] 


2.2-5 
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and 


(^(r,  0)  =  0 


8(|> 


*1(V  T)  =  -al  ~w  (V  T) 


2.2-6 


9<J) 


3r 


(1,t)  =  -Bii  4>1  (1 ,  t)  . 


Solutions  of  equations  2.2-3  and  2.2-4  which  satisfy  the  above  con- 
ditions  are: 


<f>0(r,  t)  = 


-F0  (t)  in  S 


2.2-7 


2  r  .2  ,Fo  +  1 


4>-,  (r .  t)  =  F-j  (t)  [r  In  —  -  r  - )]  +  f2(t)  In  r  +  F3(t)  2.2-8 


where 


Bii 


Fo^  1  -  Bi .  In  a 

i  ( 


F-|  (t)  = 


da  F  a. 

O  O  0 

dx  4 


F  F  ^a 

F2(t)  =  ^~2  ^  +  BT7  ^  +  BT7^  ‘  FF'  +  Fo>  "  °a()1 


and 


f3(t)  = 


Fp  F-.  2  ] 

BTT  +  ?  [1  +  bTT  (1  +  BTT^ 

i  a.  i  i 


The  solutions  2.2-7  and  2.2-8  are  incomplete  because  they  contain 
the  time-dependent  functions  aQ  and  a-j  which  are  unknown.  However,  these 
functions  may  be  found  by  substituting  equations  2.2-1  and  2.2-2  into 


2.1-4  and  collecting  like  powers  of  a.  The  first  two  of  the  resulting 
infinite  set  of  differential  equations  are 


12 


^0 

dx 


1  +  Bi .  3(f) 
l  To 


Bi.  3r  x"o 


(a„s  T) 


and 


da 


1  1  +  B  i  -j  3^1 


dx 


Bii 


[■ 


3r 


320»o 

(a  ,  t)  +  a,  - j  (a  ,  x)] 

u  3r 


with  the  associated  initial  conditions  (from  equation  2.1-5) 


ao(0)  =  1 


and 


a-j  (0)  =  0 


Equations  2.2-7  and  2.2-9  combine  to  give 


da 


1  +  Bi. 

1  o 


dx 


Bi 


i 


or  dx  =  - 


Bii 


r+  Bii  F0  da° 


which  integrates  subject  to  2.2-11  to  give 


x  = 


Bi.  a,/  In  a  ?  Bi .  +  2 

+(i-.02)4ir 


}  . 


2.2-9 


2.2-10 


2.2-11 


2.2-12 


2.2-13 


Now  using  equations  2.2-7  and  2.2-8  in  equation  2.2-10  gives 


. 

■ 
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da^ 

dx 


1  +  Bi.  da  F  2  9  F  2 

-TV1 <-£  t-f  o  +  +  -$-  d  +  fo) 


4a. 


*7  0  +  BT-  0  +  Bl-))] 


Bii 


c  2 

F  a-i  a-. 

0  1  +  1  0} 

O  T  0-1 


Rearranging 


da,  1  +  Bi.  F„  F„  da 

dir+  ai  c-rrr  (72'^)/(^)] 

o  o 


1  +  Bi,  F 


2  3 

F_  F  J 


_ 1  r_£  +  J9L  +  '  o 

Bi.  L  2  4 


4a 


VO  +  BT-(1  +  pF-))] 


Bii 


Bii 


da. 


Substituting  for  ^rom  equation  2.2-12  this  becomes 


da- 


1  -  F 


da 


+  (  aT-^)  al  - 


1  +  Bi_.  F 
Bi.  LT 


2  3 

F  F  J 

i  r_o  + 

I  'T  ^  A  A 


4a. 


^  (1  +BT-0  +  BT^))]  • 


Bii 


2.2-14 


Formal  integration  of  equation  2.2-14  subject  to  2.2-11  then  gives 


1  -  F. 


-/ 


1  +  Bi. 
-Ti— 


F^da 

e  a°  °  da 


da 


°  +  "~f~  +  T  d  +  bTT  d  +  1T7^ 


2  .3 


9Cf  -  J  i 


»  K  t 


14 


TIME  (r) 


FIG.  2.2:  ZEROTH-ORDER  SOLUTION  FOR  INTERFACE  GROWTH 
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FIG.  2.3:  FIRST  ORDER  CONTRIBUTION  TO  INTERFACE  GROWTH 
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FIG.  2.4(a):  THE  EFFECT  OF  SENSIBLE  HEAT  AND  CONVECTION 

ON  I CE  GROWTH 
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FIG.  2.4(b):  THE  EFFECT  OF  SENSIBLE  HEAT  AND  CONVECTION 

ON  ICE  GROWTH 
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FIG.  2.4(c):  THE  EFFECT  OF  SENSIBLE  HEAT  AND  CONVECTION 

ON  ICE  GROWTH 
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Equations  2.2-13  and  2.2-15  may  be  used  to  express  as  functions  of 
time  both  aQ  and  a-j ;  and,  for  a  given  Stefan  number,  the  first  order 
(i.e.  aQ  +  aa-j )  approximation  to  the  interface  location.  Figures  2.2 
and  2.3  reveal  the  time  dependency  of  aQ  and  a-j  for  various  Biot  num¬ 
bers.  Figure  2.4  reveals  the  effect  of  Stefan  number  and  Biot  number 
on  interface  growth. 


2.3  Non-Uniform  Convection 

The  results  of  section  2.2  may  be  used  to  predict  interface  shapes 
generated  by  circumferential  variations  in  the  heat  transfer  coefficient 
i.e.,  in  the  local  Biot  number.  The  results  will  have  validity  as 
long  as  a  close  approximation  to  axi-symmetric  conditions  exist  in  the 
ice  shell  (see  Appendix  A). 

Figures  2.5  and  2.6  represent  changes  in  the  interface  size  and 
shape  for  two  cases  of  asymmetric  convective  cooling.  The  Biot  number 
variations  were  obtained  from  established  data  [12]  for  transverse  flow 
over  an  isothermal  cylinder  at  various  values  of  the  external  Reynolds 
number. 

During  the  growth  of  the  ice  shell  any  water  flowing  in  the  tube 
must  experience  an  increased  resistance  to  flow,  due  to  the  reduced 
cross-sectional  area  of  the  liquid  region.  This  poses  the  question  of 
how  interface  growth  changes  the  pressure  drop. 
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REYNOLDS  NUMBER  =  23 
STEFAN  NUMBER  =  0.10 


LOCAL  NUSSELT  NUMBER 


FIG.  2.5:  ICE  GROWTH  IN  A  TRANSVERSE  AIR  H.0W:  R£  -  23 
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REYNOLDS  NUMBER  =  110,000 
STEFAN  NUMBER  =  0.10 

Bi;  =  1.48 


FIG.  2.6:  ICE  GROWTH  IN  A  TRANSVERSE  AIR  FLOW:  R£  =  11x10 


CHAPTER  III 


PRESSURE  DROP  ANALYSIS 


3.1  Governing  Equations 

The  fluid  velocity  is  governed  by  the  momentum  equation  assuming 
fully  developed  flow.  That  is 
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the  equation  may  be  written 
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3.1-1 


It  is  clear  from  figures  2.5  and  2.6  that  the  interface  does  not 
generally  possess  a  regular  shape  and  therefore  it  is  necessary  to  solve 
equation  3.1-1  within  arbitrarily  shaped  boundaries.  There  are  several 
available  techniques  [13,  14,  15]  for  such  problems.  Sparrow  and 
Haji-Sheikh  use  sophisticated  but  cumbersome  orthonormalizations  to 
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FIG.  3.1:  IMBEDDING 
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give  solutions  for  arbitrary  shapes.  Tao  uses  conformal  mapping  to 

permit  solutions  for  any  shapes  for  which  suitable  mappings  can  be 

found.  In  view  of  the  difficulties  and  limitations  involved  in  these 

and  other  techniques  it  is  proposed  to  develop  an  alternative  semi- 

* 

analytic  method  ,  similar  to  that  employed  by  Boley  [14],  which  has 
the  merit  of  relative  simplicity. 

This  method  treats  the  region  of  interest  as  an  interior  part  of 
a  larger  region  of  any  convenient  size  and  shape.  The  solution  to  equa¬ 
tion  3.1-1  is  then  found  for  the  larger  region  and  adjusted  through  the 
outer  boundary  values,  until  values  at  the  surface  of  the  inner  region 
satisfy  the  boundary  conditions  for  the  original  problem. 

3.2  Solution  by  the  Imbedding  Technique 

Consider  the  solution  of  equation  3.1-1  in  the  interior  region 
shown  in  figure  3.1.  This  region  represents  the  area  available  for 
liquid  flow  and  is  imbedded  in  the  larger  region  bounded  by  the  circle 
of  radius  R*.  This  circular  outer  boundary  is  divided  into  a  number  of 
equal  segments.  The  arc  segment  is  assumed  to  have  a  velocity  U . , 
the  magnitude  and  sign  of  which  is  to  be  determined. 

The  velocity  distribution  resulting  solely  from  the  movement  of 
one  arc  is  readily  found,  and  by  superposition  it  is  found  that  for  m 
equally  wide  segments  the  velocity  field  is  given  by 

u(r*,  0)  =  "  9x"  (1  “  r*)  “  j  Uj  ^j  ^r**  3.2-1 

3  ■ 


Analytic  in  formulation  but  numerical  in  execution. 
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where 
fj(r*,  0) 
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00 


n 

r*  sin 
irn 


cos  n  (2tt  9  -  (j  -  j)  —■) 


It  now  remains  to  choose  the  U.  such  that  the  velocity  at  the  ice- 

J  J 

water  interface  is  zero.  Thus,  for  any  co-ordinate  point  (a*^,  9^)  on 
the  interface 

.1  Uj  f.  (a*k,  0k)  =  -  If  (1  -  a*^)  .  3.2-2 

J  1 

It  is  evident  that  a  set  of  m  such  points  on  the  interface  will  yield 
a  set  of  m  linear  algebraic  equations  relating  the  arc  velocities  IK 
to  any  given  set  of  interface  radii.  The  equations  may  be  solved  read¬ 
ily  to  determine  a  set  of  U.  which  will  produce  a  zero-veloci ty  contour 

<1 

passing  through  the  chosen  co-ordinate  points. 

As  written,  equation  3.2-2  cannot  be  solved  unless  the  pressure 
gradient  is  specified.  To  circumvent  this  difficulty  we  define 

G(r*,  0)  =  -u(r*,  0)/|^ 

thus  permitting  the  calculation  of  U./ff  instead  of  IK .  The  variation  of 
may  now  be  found  from  the  volumetric  flow  rate  given  by 

dX 
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FIG.  3.2(a):  EFFECT  OF  SENSIBLE  HEAT  AND  CONVECTION  ON 

PRESSURE  GRADIENT 


PRESSURE  GRADIENT 
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FIG.  3.2(b):  EFFECT  OF  SENSIBLE  HEAT  AND  CONVECTION  ON 

PRESSURE  GRADIENT 


PRESSURE  GRADIENT 
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FIG.  3.2(c):  EFFECT  OF  SENSIBLE  HEAT  AND  CONVECTION  ON 

PRESSURE  GRADIENT 


PRESSURE  GRADIENT 
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FIG.  3.3:  PRESSURE  STEEPENING  IN  A  TRANSVERSE 
AIR  FLOW 
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For  a  constant  flow  rate  this  may  be  equated  to  the  ice-free  value 


/9P\  £ 

4y  l3X;ni  *  2  * 


Therefore 


Throughout  the  above  analysis  R*  is  always  taken  at  a  value  great 
enough  to  ensure  convergence  of  f.(r*,  0)  but  small  enough  to  permit 
motion  of  the  arcs  around  the  fictitious  boundary  to  exert  sufficient 
influence  at  the  actual  boundary.  It  was  found  worthwhile  re-establ ish- 
ing  R*  for  each  value  of  t  at  which  the  pressure  drop  was  calculated. 

Figure  3.2  reveals  the  steepening  of  the  pressure  gradient  for 
representative  values  of  the  Stefan  and  Biot  numbers  under  uniform  external 
convection.  Figure  3.3  shows  similar  curves  for  each  of  the  two  trans¬ 
verse  flows  incorporating  non-uniform  convection  under  the  same  condi-. 
tions  as  previously  considered. 
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CHAPTER  IV 


DISCUSSION  OF  THEORETICAL  RESULTS 


The  effect  of  sensible  heat  and  external  convection  on  ice  growth 
under  uniform  conditions  is  shown  in  figures  2.4  (a),  (b)  and  (c).  The 
effect  of  sensible  heat  is  revealed  through  comparison  of  the  solid 
curves,  corresponding  to  the  inclusion  of  sensible  heat  effects,  and 
the  broken  curves,  corresponding  to  the  neglect  of  sensible  heat  in 
the  ice,  that  is,  Stefan  number  equal  to  zero.  It  is  important  to  note 
that  Ste.  =  0  in  this  context  does  not  imply  that  (Tf  -  T  )  =  0,  rather 

I  I  v 

that  L  =  °°  and/or  Cp^  =  0.  Mathematically  it  represents  the  complete 
neglect  of  sensible  heat  effects;  that  is  the  quasi-steady  state  situa¬ 
tion  treated  by  London  and  Seban  [3]. 

The  range  of  Biot  numbers  for  which  the  results  are  plotted  is 
effectively  the  full  range  of  permissible  Biot  numbers.  For  very  large 
Biot  numbers  the  solution  is  asymptotic  to  the  Dirichlet  type  of  bound- 

4 

ary  condition  (Bi ^  =  °°).  Thus  the  curve  plotted  for  Bi.  =  10  is  accur¬ 
ate  for  all  very  large  Biot  numbers  (Bi.  >  10^).  In  terms  of  real  time, 
the  solutions  for  very  small  Biot  numbers  are  not  asymptotic,  but  with 


reference  to  the  normalized  time  scale  where  t  includes  Bi . ,  the  curves 

are  asymptotic  to  the  solution  for  an  adiabatic  boundary  condition 

-4  , 

(Bi .j  =  0).  Thus  the  curve  plotted  for  Bi\  =  10  is  accurate  for  all 
smaller  values  of  Bi.  as  it  tends  to  zero.  This  curve  is  given  by  the 
limit  of  equation  2.2-13  as  Bi .  approaches  zero,  noting  that  under  the 
same  limit  aQ  -*  a.  Thus 

a  =  (1  -  2t)1/2 

Bi.  0 


. 
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which  is  seen  to  be  independent  of  Ste^.  Physically  this  means  that 
the  ice  forms  so  slowly  that  the  temperature  profile  within  the  ice 
does  not  differ  from  the  steady  state  profile,  and  sensible  heat  ef¬ 
fects  are  negligible.  The  influence  of  sensible  heat  effects  is  most 
noticeable  for  large  Biot  numbers  although  it  is  evident  from  figures 
2-4  that  this  influence  is  generally  small.  Since  the  results  are  gen- 
erally  accurate  to  a  an  upper  limit  on  the  Stefan  number  of  0.3  (Fig. 
2.4(b))  ensures  that  inaccuracies  of  10%  will  not  be  exceeded  even  as 
the  Biot  number  approaches  infinity.  For  the  ice-water  system  Stefan 
numbers  greater  than  0.3  represent  coolant  temperatures  less  than 
-54°F  which,  although  possible,  occur  infrequently  in  engineering  situa¬ 
tions. 

The  minima  revealed  in  figures  2.4  have  no  physical  meaning  but 
result  from  a  singularity  in  the  first  order  solution  for  a-j .  Since 
all  higher  order  solutions  in  a  perturbation  expansion  are  dependent 
on  the  zeroth  order  solution,  it  is  to  be  expected  that  the  higher 
order  solutions  exist  only  in  the  domain  where  the  zeroth  order  solution 
exists.  In  this  case  the  zeroth  order  solution,  as  given  by  equation 
2.2-13,  exists  only  for  1  >  aQ  >  0.  Now  as  a  first  approximation 

a  =  aQ  +  a  a] 

where  both  a  and  a-j  are  positive.  Thus,  if  a  is  to  approach  zero  it  is 
apparent  that  aQ  must  become  negative.  The  solutions  for  aQ  and  a-j  do 
not  exist  for  a^  less  than  zero,  and  in  fact  a-j  becomes  infinite  as  a^ 
approaches  zero.  This  results  in  the  minima  of  figuves  2.4.  Results 
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in  the  region  of  the  minima  are  not  considered  valid. 

If  the  tube  is  immersed  in  a  transverse  flow  the  convective  bound¬ 
ary  condition  ( Bi ^ )  varies  circumferentially  in  a  manner  which  depends 
on  the  external  Reynolds  number  and  on  the  Prandtl  number  of  the  cooling 
fluid.  If  conditions  are  such  as  to  provide  a  close  approximation  to 
ax i -symmetry, i ce  growth  can  be  calculated  using  the  results  of  Chapter 
II  with  the  local  Biot  number  for  any  azimuthal  location.  Thus, an  in¬ 
terface  shape  which  approximates  the  actual  interface  shape  can  be 
generated  at  any  time.  Since  this  method  ignores  the  effects  of  azi¬ 
muthal  conductions  in  the  ice, a  good  measure  of  the  error  is  to  com¬ 
pare  the  magnitudes  of  the  radial  and  azimuthal  conduction  terms  in  the 
interface  shape  calculated  (see  Appendix  A).  Equation  A-2  in  Appendix 
A  reveals  that  the  error  approaches  zero  for  very  larges  values  of  the 
Biot  number  and  for  very  small  azimuthal  variations  as  would  be  expected. 
Figure  A.l  gives  results  for  which  the  magnitude  of  the  azimuthal  con¬ 
duction  term  is  10%  of  that  of  the  radial  conduction  term. 

Figures  2.5  and  2.6  show  ice-water  interface  profiles  for  ice 
formation  in  a  cooling  air  flow  at  two  different  Reynolds  numbers. 

Figure  2.5  describes  the  lower  Reynolds  number  situation  which  has  a 
low  mean  Biot  number.  As  is  expected,  moderate  variations  in  the  local 
Biot  number  result  in  significant  variations  in  ice  thickness.  The 
radial ly-symmetric  approximation  seems  reasonably  valid  in  this  case. 
Figure  2.6  presents  similar  results  for  a  much  higher  Reynolds  number 
and  higher  mean  Biot  number.  Although  the  variations  in  ice  thickness 
are  less  pronounced  in  this  case,  it  can  be  readily  seen  from  figure 
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2.6  that  as  lower  interface  radii  are  encountered  the  radial ly-symmet- 
ric  assumption  is  doubtful.  This  is  in  complete  agreement  with  the 
predictions  of  Appendix  A. 

One  important  engineering  consideration  in  the  formation  of  ice 
in  a  pipe  is  the  effect  of  the  reduced  area  for  liquid  flow  in  the 
pipe.  Thus,  for  a  fixed  flow  rate,  the  pressure  gradient  in  the  fluid 
will  experience  an  increase  as  the  flow  area  is  reduced.  Figures  3.2 
(a,  b,  c)  show  the  pressure  gradient  increase  for  the  cases  of  uniform 
external  convection  of  figures  2.4  (a,  b,  c).  Under  uniform  convection, 
the  pressure  gradient  is  inversely  proportional  to  the  forth  power  of 
the  interface  radius.  Figure  3.3  presents  the  results  for  the  two 
cases  of  non-uniform  convection  as  described  in  figures  2.5  and  2.6. 

It  is  apparent  that  for  the  time  interval  considered  there  is  little 
difference  between  the  pressure  gradient  calculated  for  non-uniform 
convection  using  the  interface  profiles  of  section  2.3, and  that  calcul¬ 
ated  assuming  uniform  convection  at  the  mean  Biot  number  of  the  trans¬ 
verse  flow. 

The  analysis  assumes  that  the  water  is  not  superheated.  However 
the  results  may  be  considered  valid  for  situations  where 
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Under  conditions  where  the  water  is  only  heated  by  viscous  dissipation 
the  temperature  in  the  water  is  given  approximately  by 
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from  which  it  follows  that  viscous  dissipation  may  be  ignored  if 
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that  is 

uav2  «  L<!r>  £t>  stei  Bii  • 
w  pi 

Using  the  appropriate  properties  for  ice  and  water  it  can  be  shown  that 
for  a  pipe  of  radius  A  (in  feet),  the  effect  of  viscous  dissipation  may 
be  neglected  if 

Re  <  107  A  (Ste.  Bi.)1/2  . 


It  is  clear  that  dissipation  will  not  normally  be  important  unless  Ste. 
Bi.,  or  a,  is  extremely  small  and/or  the  water  velocity  is  very  high. 
Physically  a  very  small  value  of  a  implies  either  a  very  small  tempera¬ 
ture  difference  or  a  very  small  convective  heat  transfer  coefficient. 
Either  of  these  situations  suggests  very  slow  ice  growth  in  which  case 
even  small  energy  contributions  from  viscous  dissipation  will  have  a 
significant  effect.  For  a  fixed  flow  rate  very  high  water  velocities 
will  undou'otably  occur  as  the  ice-water  interface  approaches  the  tube 
axis.  Thus  results  can  be  inaccurate  when  the  water  occupies  only  a 
small  region,  e.g.  near  the  tube  centerline.  If  the  pressure  drop  is 
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constant,  dissipation  will  be  less  important  unless  it  cannot  be  initi¬ 
ally  ignored. 

The  analysis  ignores  circumferential  conduction  in  the  tube  wall 
since  it  is  reasonable  to  assume  w  «  R-|  and  isotherms  are  expected  to 
be  roughly  circular.  Therefore  azimuthal  temperature  gradients  will 
usually  be  small  and  will  not  generate  significant  heat  fluxes  in  the 
azimuthal  direction.  In  the  radial  direction 


k  21 1  =  k  —  I 

i  3R  1  R1  m  8R  1  R] 

or  approximately 

AT.  AT 

k  . .  1  ~  k  _ ~ 

Ki  R]  "  m  w 

and  hence 

A^"m  _  ki  w 

"  km  Ri  ' 

Thus,  assuming  the  tube  is  metallic  it  follows  that  AT^  «  AT.  and  does 

not  significantly  alter  the  assumed  boundary  conditions.  Since  c  <  c 

the  effect  of  sensible  heat  in  the  metal  is  likewise  insignificant. 

Although  the  system  of  ice  and  water  has  been  central  throughout 

the  analysis  the  discussion  is  not  limited  to  such  a  system.  The  only 

restriction  to  be  placed  on  any  physical  situation  is  that  a.  must  be 

2 

less  than  unity,  noting  that  the  solution  is  in  error  by  a  . 
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PART  II 


EXPERIMENTAL  INVESTIGATION 


CHAPTER  V 


EXPERIMENTAL  APPARATUS 

The  experimental  apparatus  employed  was  essentially  that  used  by 
Freeborn  [11]  with  certain  modifications.  It  was  considered  as  two 
systems;  one  for  controlling  the  external  boundary  conditions  and  the 
other  for  the  internal  fluid  flow  (see  fig.  5.1), 

The  external  boundary  condition  was  controlled  by  a  closed  system 
primarily  composed  of  a  reservoir,  pump  and  piping  which  supplied  the 
test  section  cooling  jacket  with  a  constant  flow  of  the  cooling  agent 
(equal  quantities  of  water  and  ethylene  glycol).  The  coolant  was  stored 
in  a  3.5  cu.  ft.  reservoir  equipped  with  a  cooling  coil  connected  to  a 
vapor  compression  refrigerati on  unit.  The  refrigeration  capacity  of 
one  ton  enabled  the  reservoir  to  be  maintained  at  a  constant  temperature 
of  -15°F.  The  reservoir  supplied  the  cooling  jacket  through  a  one-inch 
diameter  copper  tubing  recirculating  system,  containing  a  motorized 
three-way  valve  and  pump. 

A  resistance  bulb  thermometer,  (L  7033A  Balco)  located  about  one 
foot  downstream  of  the  three-way  valve,  was  used  with  a  resistance  therm¬ 
ometer  controller  (Honeywell  R  70870)  to  actuate  the  three-way  valve 
(Honeywell  Series  1616  three-way  valve  with  M904  Modutol  Motor).  The 
pump  recirculated  the  coolant  from  the  jacket  to  the  main  reservoir 
and  to  the  three-way  valve  which  mixed  some  returning  coolant  with 
-1 5°f  coolant  direct  from  the  reservoir  to  maintain  the  temperature  set 
on  the  resistance  thermometer  controller.  The  reservoir  and  piping 
were  insulated  to  prevent  heat  leakage  to  the  system. 
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1  THREE-WAY  VALVE 

3  GATE  VALVE 


FIG.  5.1:  SCHEMATIC  OF  EXPERIMENTAL  APPARATUS 
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The  water  flow  was  governed  by  an  open  system  consisting  of  two 
constant-head  storage  tanks,  supply  and  discharge  piping,  motorized 
three-way  valve,  flow  meter,  entry  length  and  test  section. 

Water  was  stored  in  the  larger  header  tank  (9  cubic  feet)  where 
it  was  cooled  by  a  coil  around  which  an  icebank  of  fixed  thickness  was 
allowed  to  form.  Thus  the  water  was  maintained  at  a  constant  tempera¬ 
ture  near  the  freezing  point.  The  coil  was  connected  to  a  vapor-com¬ 
pression  refrigeration  unit  with  a  one-half  ton  capacity.  The  water 
in  the  second  and  smaller  header  tank  (1  cubic  foot)  was  maintained  at 
the  water  main  temperature.  Water  from  both  tanks  was  piped  to  the 
motorized  three-way  valve  where  the  "cold"  water  was  trimmed  to  the 
desired  temperature  by  mixing  with  the  "warm"  water  from  the  small  tank. 
The  water  was  fed  through  one  inch  diameter  copper  tubing,  through  a 
flow  meter  (Fisher  and  Porter  Ratosight)  and  controlling  gate  valve 
into  a  six  foot  section  of  three  inch  diameter  copper  tubing;  the  last 
three  feet  of  which  constituted  the  test  section.  This  section  was 
surrounded  by  the  cooling  jacket  and  the  complete  system  was  thermally 
insulated. 

The  complete  test  section  and  inlet  piping  were  mounted  on  a  frame 
which  enabled  the  test  section  to  be  rotated  to,  and  held  at,  any  in¬ 
clination.  To  allow  freedom  of  operation  the  inlet  water  connection 
was  made  through  a  swivel  joint  and  the  coolant  system  incorporated 
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CHAPTER  VI 


INSTRUMENTATION 

The  prime  instrumentation  problem  involved  measuring  the  ice  thick¬ 
ness  inside  a  long  vertical  pipe  completely  surrounded  by  an  insulated 
cooling  jacket.  In  the  general  case  the  ice  thickness  has  both  axial 
and  azimuthal  variation. 

Gort  [9]  measured  the  ice  thickness  with  a  stainless  steel  probe 
from  which  mechanical  arms  were  deployed  to  contact  the  ice  water  in¬ 
terface.  The  device  was  only  suitable  for  radially  symmetric  situations 
although  it  did  permit  an  axial  traverse. 

Freeborn  [11]  was  successful  at  measuring  ice  thickness  with  an 
ultrasonic  technique.  However,  attempts  at  obtaining  axial  profiles 
met  with  little  success  except  in  regions  where  axial  variations  in 
the  ice  thickness  were  very  small.  The  ultrasonic  technique  has  the 
very  definite  advantage  of  measuring  ice  thickness  directly  as  opposed 
to  measuring  the  distance  from  the  probe  to  the  ice-water  interface. 

The  latter  requires  knowledge  of  the  exact  location  of  the  probe  rela¬ 
tive  to  the  pipe  centerline. 

In  order  to  study  both  axial  and  azimuthal  variations  in  the  ice 
thickness  a  new  finger-type  mechanical  probe  was  designed  (fig.  6.1). 

The  device  was  intended  to  operate  as  a  continuous  reading  profilometer 
to  be  drawn-  axially  through  the  pipe.  Variations  in  the  ice  thickness 
cause  radial  movement  of  the  sensing  arm  which  results  in  a  correspond¬ 
ing  axial  movement  of  the  transducer  core.  This  causes  a  change  in  the 
electrical  signal  from  the  voltage  transducer  (Ametek  400-3K-9EL)  pro¬ 
portional  to  the  radial  displacement.  Calibration  of  the  device  with  a 
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FIG.  6.1:  SCHEMATIC  OF  FINGER  PROF  I  LOME  TER 
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micrometer  screw  revealed  an  effective  resolution  better  than  0.005 
inches  and  readings  were  repeatable  to  within  this  limit. 

The  axial  position  of  the  probe  was  controlled  by  a  reversing  D.C. 
motor  driving  a  lead  screw  to  which  the  one-half  inch  diameter  probe 
was  attached.  The  position  of  the  probe  was  then  monitored  by  a  forty- 
turn  potentiometer:  a  change  in  potential  from  a  prescribed  datum  in¬ 
dicated  the  distance  of  the  probe  from  the  test  section  entry.  The 
signals  from  the  potentiometer  and  the  voltage  displacement  transducer 
were  then  used  to  drive  an  X-Y  recorder  with  the  ice  thickness  plotted 
on  the  vertical  axis  and  the  axial  position  on  the  horizontal  axis. 

By  this  means  an  axial  profile  of  the  ice  thickness,  for  any  azimuthal 
location  at  a  particular  time,  could  be  generated. 

To  generate  azimuthal  profiles  the  probe  was  set  at  the  desired 
axial  locations  and  with  a  time  base  on  the  X-Y  plotter  the  probe  was 
rotated  through  180  degrees  with  a  five  second  pause  at  each  azimuthal 
location. 

In  order  to  maintain  the  desired  coolant  temperature,  a  resistance 
bulb  thermometer  was  located  in  the  coolant  recirculating  line,  as 
previously  mentioned.  The  signal  from  this  thermometer  was  transmitted 
to  the  temperature  control  mechanism  where  it  was  compared  to  a  pre¬ 
set  control  value.  The  error  signal,  if  any,  energized  the  motor  con¬ 
trolling  the  three-way  valve  movement,  thus  positioning  the  valve  for 
a  higher  proportion  of  either  the  -15°F  coolant  or  the  warmer  recircu¬ 
lating  coolant.  The  three-way  valve  operation  for  control  of  the  water 
temperature  was  similar,  except  that  a  copper-constantan  thermocouple 
sensor  was  used  and  the  valve  movement  was  controlled  manually  by  use 
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of  a  three-way  toggle  switch. 

Further  instrumentation  consisted  of  copper  constantan  thermo¬ 
couples  situated  in  the  water  inlet  region,  six  and  nine  inches  before 
the  test  section  entry,  and  in  the  coolant  recirculation  piping  immed¬ 
iately  downstream  from  the  cooling  jacket.  In  an  attempt  to  measure 
the  temperature  gradient  across  the  pipe  wall  in  the  test  section, 
thermocouples  were  located  on  the  inner  and  outer  surfaces  of  the  test 
section  in  three  locations.  Readings  were  obtained  using  either  a  milli¬ 
volt  potentiometer  (Leeds  &  Northrup  8686)  or  a  digital  voltmeter 
(Hewlett-Packard  3440A). 


CHAPTER  VII 


EXPERIMENTS 


7.1  No-Flow  Test 

Initially  a  test  was  conducted  at  zero  flow  rate  and  zero  superheat 
(a  =  0)  in  order  to  make  a  comparison  with  theory  and  other  experiments 
(see  Freeborn  [in). 

Before  the  test,  the  water  in  the  large  header  tank  was  cooled 
over  a  prolonged  period.  The  test  section  was  completely  drained  and 
set  in  the  vertical  position.  A  traverse  was  made  with  the  profi lometer, 
giving  a  datum  line  on  the  X-Y  recorder  for  the  no-ice  condition.  The 
cooling  system  recirculating  pump  was  then  started  with  the  coolant 
temperature  control  set  at  approximately  -5°F.  After  waiting  to  ensure 
the  stabilization  of  the  coolant  temperature  the  valve  to  the  test  sec¬ 
tion  was  opened  until  water  began  to  emerge  from  the  discharge  line. 

The  valve  was  then  closed  and  the  timer  started.  Axial  traverses  were 
made  at  approximately  three  minute  intervals  and  thermocouple  readings 
were  recorded  as  the  traverse  was  being  made.  At  the  conclusion  of  the 
test  the  coolant  temperature  control  was  set  above  the  freezing  point 
and  warm  water  passed  through  the  test  section,  melting  the  ice.  After 
the  ice  was  melted  the  water  was  shut  off  and  the  test  section  drained. 

7.2  Inclined  Laminar  Flow  Tests 

Next  a  series  of  tests  was  conducted  to  determine  the  influence  of 
flow,  superheat  and  angle  of  inclination  on  the  rate  of  ice  formation. 
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With  the  coolant  temperature  control  set  above  the  freezing  point, 
the  flow  was  directed  through  the  test  section  and  adjusted  to  get  the 
desired  flow  rate  and  superheat  ratio  with  the  test  section  set  at  the 
desired  inclination.  The  flow  was  then  discharged  before  reaching  the 
test  section  and  the  test  section  was  prepared  as  described  in  the  no¬ 
flow  test.  With  the  system  in  a  steady  condition  the  flow  was  switched 
back  through  the  test  section  and  the  timer  started.  Axial  profiles 
and  temperatures  were  recorded  as  before.  In  addition,  some  azimuthal 
traverses  were  made  with  the  probe  at  a  fixed  axial  location.  Periodi¬ 
cally,  the  flow  rate  was  checked  by  measuring  the  amount  of  water  dis¬ 
charged  over  a  timed  interval.  At  the  end  of  the  test  the  ice  was 
melted  and  the  section  drained  as  previously  described. 
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CHAPTER  VIII 


DISCUSSION  OF  EXPERIMENTAL  RESULTS 

8.1  No-Flow  Test 

This  experiment  was  to  be  conducted  without  flow  and  with  zero 
superheat  as  previously  mentioned.  Ideally,  the  situation  is  that  of 
a  liquid  flowing  at  its  freezing  temperature  in  a  pipe  cooled  by  uni¬ 
form  external  convection  (i.e.  fully  developed  ice  formation). 

In  actual  fact  the  water  temperature  was  approximately  35°F  with 
a  superheat  ratio  of  0.09.  Figure  8.1  shows  the  results  of  the  axial 
traverses  for  this  test.  There  are  several  reasons  for  the  axial  vari 
ations  in  the  ice  thickness.  Near  the  entry  region  the  effects  of 
axial  conduction  in  the  ice  would  tend  to  reduce  the  rate  of  ice  forma 
tion  as  can  be  seen  from  the  graph  near  z  =  0.  Secondly,  since  the 
water  is  not  initially  at  the  freezing  temperature  there  must  be  some 
effect  due  to  natural  convection  and  this  effect  will  have  an  axial 
dependency  contributing  to  axial  variations  in  the  ice  thickness  at 
any  time.  A  third  contributing  factor  is  the  probability  that  the  ex¬ 
ternal  heat  transfer  coefficient  is  non-uniform  in  the  axial  direction 

It  should  be  noted  that  in  figure  8.1  the  graph  scaling  is  such 
that  ice  thickness  variations  are  greatly  exaggerated.  Figure  B. 1  of 
Appendix  B  shows  the  same  curve  in  the  proper  perspective.  Here  it 
becomes  apparent  that  the  experiment  is  in  fact  very  close  to  fully 
developed  ice  formation. 

Figure  8.2  shows  the  test  results  at  the  axial  location  where 
Freeborn  [11]  obtained  results  with  the  ultrasonic  profilometer.  Also 
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FIG,  8.1:  AXIAL  PROFILES 
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FIG.  8.2:  COMPARISON  WITH  PREVIOUS  RESULTS 
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shown  are  the  results  of  Freeborn's  no-flow,  zero  superheat  test  and 
the  corresponding  results  from  Gort's  [9]  experiments.  There  is  very 
good  agreement  between  the  results  of  all  three  methods. 

The  solid  and  broken  curves  in  figure  8.2  are  the  zeroth  order 
and  first  order  correct  theoretical  solutions  for  fully  developed  ice 
formation  as  developed  earlier  in  this  thesis.  The  first  order  solu¬ 
tion  is  terminated  short  of  the  extremum  and  is  seen  to  be  in  excellent 
agreement  with  all  of  the  experimental  results. 

An  attempt  was  made  to  determine  the  external  heat  transfer  coef¬ 
ficient  by  measuring  the  temperature  gradient  across  the  pipe  wall. 
There  were  inner  and  outer  thermocouples  installed  at  three  locations 
but  complete  failures  occurred  at  two  of  the  locations.  At  the  third 
location  there  were  two  adjacent  thermocouples  on  the  inside  wall  of 
the  pipe  but  their  readings  were  in  considerable  disagreement.  Con¬ 
sequently  this  method  of  establishing  the  heat  transfer  coefficient  was 
abandoned.  Because  of  the  excellent  agreement  between  the  experimental 
results  of  this  test  and  the  corresponding  results  of  Freeborn  [11], 
the  external  heat  transfer  coefficient  was  assumed  to  be  the  same  as 
in  Freeborn's  work. 

8.2  Inclined  Laminar  Flow  Tests 

Tests  were  run  at  four  different  test  section  inclinations  for 
laminar  flows  with  superheat.  The  results  of  the  tests  are  presented 
in  figures  8.3  through  8.13.  In  each  case  the  axial  variables  are 
nondimensional ized  to  give 
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z  =  Z/Pr  Re 

The  inclination  of  the  test  section  is  y  given  in  degrees  from  the  hori¬ 
zontal.  In  each  inclined  test  the  axial  profiles  were  taken  on  the 
lower  wall  of  the  pipe. 

Initially,  it  had  been  planned  to  keep  the  Reynolds  number  and 
superheat  ratio  constant  throughout  the  tests  and  explore,  through  the 
angle  of  inclination,  the  effects  of  natural  convection.  Great  diffi¬ 
culty  was  experienced  in  controlling  the  superheat  ratio  due  to  the 
time  lag  between  operation  of  the  three-way  mixing  valve  and  correspond¬ 
ing  response  in  the  inlet  water  temperature.  Similarly,  some  difficulty 
was  experienced  in  maintaining  the  desired  Reynolds  number  in  each  test. 
The  resulting  scatter  in  these  two  parameters  makes  a  strict  analysis 
of  the  results,  with  respect  to  the  effect  of  inclination,  very  diffi¬ 
cult.  However  the  results  allow  some  general  conclusions  to  be  drawn 
about  the  nature  of  ice  formation  under  these  conditions. 

Figures  8.3,  8.5,  8.8  and  8.11  are  the  axial  profiles  at  various 
times  for  each  test.  In  each  case  they  reveal  the  existence  of  an 
asymptotic  steady  state  condition  near  the  entry  region.  Further  down¬ 
stream  where  the  heat  had  been  largely  removed  from  the  water,  the  pro¬ 
files  indicate  that  the  situation  was  still  a  transient  one.  The  ice 
thickness  is  expected  to  increase  monotonical ly  with  increasing  z  for 
the  vertical  case  (y  =  90°),  and  although  the  inclined  cases  are  more 
difficult  to  analyze,  a  similar  variation  may  be  expected  in  these 
cases  too.  However  the  axial  profiles  for  small  times  reveal  that  the 
ice  thickness  is  not  monotone.  The  strong  similarity  between  the  first 
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profiles  in  all  four  tests  suggests  axial  variations  in  the  heat  trans¬ 
fer  coefficient  and/or  the  mean  coolant  temperature.  This  is  supported 
by  the  fact  that  the  exit  region  for  the  coolant  in  the  cooling  jacket 
was  located  near  the  water  entry  region  on  the  same  side  of  the  test 
section  as  the  profiles  were  obtained.  This  suggests  that  improved 
circulation  in  the  cooling  jacket  near  the  exit  results  in  a  higher 
Biot  number  which  accounts  for  a  thickening  in  the  ice  profile  in  that 
region.  At  later  times  the  effect  of  axial  variations  in  the  Biot  num¬ 
ber  is  reduced  due  to  axial  conduction  in  the  ice  layer,  and  the  pro¬ 
files  reveal  the  expected  monotonicity. 

Figures  8.4,  8.6,  8.9  and  8.12  are  simply  cross-plots  from  the 
corresponding  axial  profiles  at  various  axial  locations.  Once  again  the 
asymptotic  steady  state  situation  is  revealed  for  small  values  of  z  and 
the  purely  transient  nature  is  exhibited  for  larger  values  of  z.  The 
ice  thickness  at  any  axial  location  is  expected  to  be  monotonical ly  in¬ 
creasing  with  time.  In  each  of  the  tests  it  is  observed  that  this  is 
not  the  case  for  the  smallest  value  of  z  plotted.  This  is  due  to  the 
fact  that  when  water  first  enters  the  test  section  the  entire  section 
has  been  cooled  to  the  temperature  of  the  coolant.  Initially  therefore 
the  boundary  condition  is  a  Dirichlet  type  and  the  ice  forms  more  rapid¬ 
ly  than  it  would  under  the  convective  boundary  condition.  For  small  z 
this  initial  growth  exceeds  the  eventual  steady-state  thickness  and, 
as  this  quantity  of  ice  cannot  be  supported  by  the  convective  boundary 
condition,  the  excess  melts  away  at  later  times.  At  larger  values  of  z, 
although  the  initial  ice  growth  is  accelerated,  the  effect  is  not  so 
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FIG.  8.4:  INTERFACE  HISTORY 
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FIG.  8.5:  AXIAL  PROFILES 
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FIG.  3.6:  INTERFACE  HISTORY 
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FIG.  8.10:  AZIMUTHAL  PROFILES  y  =  10° 
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obvious  due  to  the  fully  transient  behavior  at  these  locations  and  the 
absence  of  significant  superheat  to  cause  melting. 

Assuming  a  radially  symmetric  heat  transfer  coefficient  around  the 
test  section,  any  asymmetry  in  the  azimuthal  profiles  for  the  inclined 
tests  would  be  due  to  natural  convection  in  the  water.  At  any  particu¬ 
lar  inclination  the  degree  of  asymmetry  indicates  the  importance  of 
free  convection  for  the  given  superheat  ratio  and  Reynolds  number. 

From  the  azimuthal  profiles  of  figures  8.7,  8.10  and  8.13  it  is  appar¬ 
ent  that  free  convection  has  a  considerable  effect  under  the  experimental 
conditions.  In  each  case  the  ice  is  thicker  on  the  lower  wall  of  the 
test  section  than  on  the  top  or  sides.  With  a  constant  Reynolds  number 
and  superheat  ratio  it  is  expected  that  the  asymmetry  would  be  greatest 

for  the  horizontal  case.  Figures  8.7,  8.10  and  8.13  reveal  that  while 

the  asymmetry  for  y  =  0°  and  y  =  10°  is  considerably  greater  than  that 
for  y  =  30°,  there  is  little  difference  in  the  degree  of  asymmetry  dis¬ 
played  for  y  =  0°  and  y  =  10°.  This  is  probably  due  to  the  fact  that 

the  superheat  ratio  for  the  test  at  y  =  0°  is  significantly  less  than 

that  for  y  =  10°. 
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PART  III 


CONCLUSIONS 


CHAPTER  IX 


CONCLUDING  REMARKS 


9.1  Conclusions 

This  thesis  has  presented  a  theoretical  analysis  of  ice  formation 
in  a  convecti vely-cooled  circular  tube  in  the  absence  of  superheat  in 
the  water.  It  has  also  presented  experimental  results  for  ice  forma¬ 
tion  within  a  convecti vely-cooled  inclined  tube  in  the  presence  of  body 
forces. 

The  theoretical  analysis  accommodates  the  effect  of  sensible  heat 
in  the  ice  through  a  perturbation  expansion  in  a  typically  small  para¬ 
meter.  Generally,  sensible  heat  effects  were  found  to  be  small  and 
were  shown  to  be  of  greatest  importance  in  situations  when  the  cooling 
is  very  rapid,  i.e.  Bi ^  »  1  and/or  (T^  -  Tc)  »  1. 

The  analysis  was  applied  to  a  particular  asymmetric  convective 
cooling  situation  (a  pipe  in  a  cross-flow)  which  could  be  treated  as 
quas i -radi al ly  symmetric.  A  semi -analytic  technique  was  developed  to 
analyze  the  flow  of  water  within  the  asymmetric  ice  shells  which  de¬ 
velop  in  such  situations. 

Experiments  were  carried  out  for  ice  formation  within  a  pipe  under 
both  symmetric  (vertical)  and  asymmetric  (inclined)  test  conditions. 

A  specially-designed  finger  profilometer  enabled  both  axial  and  axi- 
muthal  profiles  of  the  ice-water  interface  to  be  ootained.  Profiles 
were  first  obtained  for  the  symmetric  case  of  ice  formation  with  no 
flow  and  no  superheat  in  the  water .  The  results  were  found  to  be  in 
excellent  agreement  with  previous  results  which  had  been  obtained  by 


. 

. 
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two  different  techniques.  From  this  agreement  it  can  be  concluded  that 
the  experiments  are  reproducible. 

Results  were  then  obtained  for  the  case  of  superheated  water  flow¬ 
ing  in  an  inclined  tube.  Ice  growth  under  these  conditions  is  not  sym¬ 
metric  and  therefore  both  axial  and  azimuthal  traverses  were  performed. 
Although  theoretical  predictions  for  ice  growth  under  these  conditions 
are  not  available  for  comparison,  the  results  appeared  consistent  with 
expectations.  For  the  region  near  the  entry  to  the  test  section  an  asy¬ 
mptotic  steady-state  interface  was  revealed  in  all  tests  in  which  super¬ 
heat  was  present.  Further  downstream  the  situation  was  found  to  be 
fully  transient  and  approaching  fully-developed  ice  formation.  The 
azimuthal  profiles  of  the  ice-water  interface  revealed  that  natural 
convection  in  the  water  is  of  considerable  importance  under  the  test 
conditions  employed.  It  was  observed  that  the  asymmetry  in  the  profiles 
was  greatest  for  inclinations  very  near  the  horizontal.  In  all  inclined 
tests  the  ice  formed  most  rapidly  on  the  lower  wall  of  the  test  section. 

9.2  Suggestions  for  Further  Work 

The  finger  profilometer  proved  very  successful  for  obtaining  both 
axial  and  azimuthal  interface  profiles.  It  is  felt  that  this  method 
could  be  used  to  obtain  results  for  a  very  extensive  experimental  study 
of  ice  formation  in  pipes.  The  existing  apparatus  would  need  improve¬ 
ments  to  allow  better  control  over  inlet  and  boundary  conditions.  For 
a  very  large  test  program  it  would  also  be  desirable  to  automate  the 
data  recording  through  the  use  of  magnetic  tape  recording  devices  com¬ 
patible  with  the  University  computing  facilities. 
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The  heat  flux  vector  for  two  dimensional  conduction  is  given  by 


q  =  -k  V*  =  qr  er  +  q0  e0 


Therefore  the  ratio  of  qA  to  q  will  be  a  measure  of  the  error  in  the 
quasi -radial ly  symmetric  analysis.  Using  the  zeroth  order  approximation 
we  find  from  equation  2.2-7 


4> 


-Bi. 

_ _ i 

1  -  Bi.  In  a 


Thus  if  Bi\  is  a  function  of  0 


_  1  3(f)  _  1  dBi  9(f)  1  da  3(f) 

%  "  r  30  r  d0  3 Bi  r  d0  3a 


Evaluating  at  r  =  a  this  gives 


I  da  1  Bi 

qQ I  a  d0  2  1  -  Bi  In  a 

a 


At  the  same  location 

i  3<f)  |  zl  (  Bi  1 
qr I a  3r'a  a  M  -  Bi  In  a' 

The  ratio  of  qQ  to  q  at  a  is  therefore  given  by 
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llii 

_  9  96 1  a  _  1  da 

04>i  ¥  dF 

3r 1  a 


The  interface  location  a  is  found  from  equation  2.2-13  which  can  be 
rearranged  to  give 


Bi  +  2 
Bi 
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1  +  Bi . 
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or  f (Bi )  =  g(a). 

Differentiating  with  respect  to  0  we  find 


d B i  /3 f  \  _  da  /  9  q  \ 
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from  which 
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Substituting  into  A-l  for  ^-yields  the  result 


E  = 


1  -  a  ,  1  ^  dBi 

9  V  9  l  rlA 

2a^  Bi.  -  Bi/  In  a  0 


A-2 


From  this  expression  it  is  possible  for  any  given  boundary  conditions 
dBi 

(Bi ,  ,  to  determine  at  what  time  (corresponding  to  a  value  of  a) 

a  given  error  will  be  exceeded.  Figure  A.l  gives  the  results  for  a 

dBi 

flux  ratio  of  E  =  0.1.  For  given  values  of  Bi  and  a  minimum  value 
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FIG.  A.l:  CURVES  FOR  EOF  TEN  PERCENT 


of  a  can  be  found  such  that  E  is  less  than  0.1. 
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APPENDIX  B 

THE  FULLY  DEVELOPED  ICE  FORMATION  PROFILE 

Figure  B.l  is  a  plot  of  the  results  of  the  no-flow,  no-superheat 
tests  as  given  in  Figure  8.1.  The  axes  correspond  to  the  actual  physi¬ 
cal  dimensions  of  the  test  section  and  the  ice  thickness. 
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FIG.  B.l:  FULLY  DEVELOPED  ICE  FORMATION 


